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(57) Abstract: The invention concerns a method for extracting a geo- 
logical horizon and related properties of a seismic representation, com- 
prising a step (100) which consists in digital modelling with continu- 
ous local seismic traces, calculating the optimal offset and defining a 
conditional neighbourhood of a reference central continuous local seis- 
mic trace; a step (101) which consists in defining a two-dimensional 
matrix whereof the line and column indices correspond to the coordi- 
nates of the geophones; a third step (102) which consists in selecting a 
seed point; a fourth step (103) which consists in determining the point 
vertically closest to the seed point and a fifth step (104) which consists 
in assigning to the point P(p,q,t) the value P(p,q,t + h ij(PqJc ), where h ijtPqX 
is the optimal offset of the neighbouring point P(i j,k), so as to estimate 
the related properties of the conditional neighbourhood thereby filling 
the two-dimensional extraction matrix of step (101). 

(57) Abrege : Un proc6d£ pour extraire un horizon g6ologique et 
des propri£t6s associees d'une representation sismique, comprend 
une 6tape 100 de modeUisation numerique par traces sismiques 
locales continues, de calcul du d£calage optimal et de definition 
du voisinage conditionnel d'une trace sismique locale continue 
"centrale" de reference ; une 6tape 101 de definition d'une 
matrice bidimensionnelle dont les indices de lignes et de colonnes 
correspondent aux coordonn6es des geophones; une troisieme 
6tape 102 de choix d'un point «graine», une quatrieme 6tape 103 
de determination du point verticalement le plus proche d'un point 
«graine» et une cinquieme 6tape 104 ou Ton attribue au point P(p,q,t) 
la valeur P(p,q,t + h ijtP q tk ), ou h ijtWfk est le decalage optimal du point 
voisin P(ij,k), de maniere a estimer les propri6t£s associees du 
voisinage conditionnel et ainsi remplir la matrice bidimensionnelle 
d'extraction de T6tape 101. 
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Dispositif et produit-programme pour extraire un horizon geologiqu e et des 
propriet6s associ6es . 

^invention est relative k un procede pour extraire un horizon geologique et des 
propriety associees de leur image d6duite de donn6es sismiques. 

5 L'invention est egalement relative k un dispositif pour extraire un horizon 

geologique et des propriet6s associees de leur image deduite de donn£es sismiques. 

L'invention est enfin relative a un produit-programme d'ordinateur pour extraire 
un horizon geologique et des propri6tes associees de leur image deduite de donnees 
sismiques. 

10 On connait des proced6s d'extraction dliorizon et de leurs attributs sismiques. 

Dans ces proc6des connus, on construit une matrice d'attribut sismique 
tridimensionnelle de meme taille et de meme echantillonnage que la matrice sismique 
tridimensionnelle initiale. 

Pour attribuer k chaque point de l'horizon une valeur d f un attribut sismique, on 
15 peut proceder soit par interpolation des noeuds environnants de la matrice d'attribut 
sismique tridimensionnelle ou choisir la valeur stockee dans le nceud le plus proche de 
la matrice d'attribut sismique tridimensionnelle. 

En cas ^interpolation, le pre-traitement correspondant n^cessite un temps de 
traitement important et une quantite de memoire disponible tr&s importante ; c'est 
20 pourquoi il est g6neralement necessaire d'effectuer ce pre-traitement sur un ordinateur 
puissant. 

Dans le cas du choix de l'attribut du noeud le plus proche, ce pre-traitement n'est 
pas necessaire, et la determination de Tattribut ne peut Stre effectu6e qu'k des points 
voisins des points du reseau de la matrice sismique tridimensionnelle originate. 

25 Un inconvenient majeur de Tart anterieur est d f introduce des defauts ou 

artefacts qui induisent des erreurs verticales sur les horizons et les attributs sismiques 
et conduisent k une mauvaise interpretation geologique des mesures sismiques. 

Un premier but de Tinvention est d f ameiiorer la technique connue en proposant 
un calcul direct et plus exact, en minimisant les risques d f erreurs possibles sur les 
30 attributs en tous points de lliorizon. 
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Un deuxieme but de l'invention est de supprimer les erreurs verticales 
susceptibles de conduire h une mauvaise interpretation g6ologique. 

L'invention a pour objet un proced^ pour extraire un horizon g^ologique et des 
proprietes associees d'une representation sismique, dans lequel on construit une 
5 fonction continue Sij,k(t) par interpolation ou approximation des traces sismiques 
discretes d'une matrice sismique multidimensionnelle, ladite fonction etant designee 
comme "trace sismique locale continue", comportant les etapes suivantes : 

a), utiliser comme decalage optimal de deux traces sismiques locales continues 
voisines, la valeur de decalage rendant maximale leur fonction de correlation ; 

10 b). retenir comme voisinage conditionnel d'une trace sismique locale continue 

"centrale" Sij,k(t) le sous-voisinage consistant en des traces voisines S pq ^(t) 
correspondant a des decalages optimaux associes h des correlations Rij >pq ,k(h) 
superieures a un seuil predetermine compris entre 0 et 1 ; 

c) . construire une matrice bidimensionnelle d'extraction destinee k etre remplie 
15 par des valeurs extraites ; 

d) . choisir un point «graine» P(i j,t) et determiner le point P(i j,k) verticalement 
le plus proche ; 

e) . estimer les proprietes associees du voisinage conditionnel et remplir la 
matrice bidimensionnelle d'extraction par des proprietes decaiees par translation de la 

20 variable courante (t) de la valeur du decalage optimal (hij, pq> k) correspondant au point 
P(i j,k) verticalement le plus proche. 

Selon d'autres caracteristiques de Tinvention : 

- on remplace successivement le contenu de la matrice bidimensionnelle de 
retape c) au cours des iterations successives d'extraction, 

25 - h chaque iteration, on utilise comme nouveaux points «graines» de retape d) 

tous les points correspondent a des emplacements de la matrice bidimensionnelle de 
retape c) remplis h Titeration precedente, 

— la propriete de subsurface extraite est Tamplitude sismique refl6chie et captee 
par des geophones, 
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- on visualise sur un 6cran de visualisation des surfaces horizon peintes avec 
des couleurs correspondant a un codage des proprietes extraites. 

L'invention est egalement relative k un dispositif pour la mise en oeuvre du 
procede selon Tinvention, comportant des moyens pour utiliser comme decalage 

5 optimal de deux traces sismiques locales continues voisines, la valeur de d6calage 
rendant maximale leur fonction de correlation, des moyens pour retenir comme 
voisinage conditionnel d'une trace sismique locale continue "centrale" S^t) le sous- 
voisinage consistant en des traces voisines S^kCt) correspondant a des d6calages 
optimaux associes a des correlations Rij >pq>k (h) superieures h un seuil predetermine 

10 compris entre 0 et 1, des moyens pour construire une matrice bidimensionnelle 
d'extraction destin6e k gtre remplie par des valeurs extraites, des moyens pour choisir 
un point «graine» P(i j,t) et determiner le point P(iJ,k) verticalement le plus proche et 
des moyens pour estimer les proprietes associees du voisinage conditionnel et remplir 
la matrice bidimensionnelle d'extraction par des proprietes decaiees par translation de 

15 la variable courante (t) de la valeur du decalage optimal correspondant au point P(i,j,k) 
verticalement le plus proche. 

Selon d'autres caracteristiques de Tinvention : 

- le dispositif comporte des moyens de memorisation et des moyens de 
visualisation de parametres sismiques determines a Taide du procede selon Tinvention. 

20 L'invention a enfin pour objet un produit-programme d'ordinateur, comportant 

des elements de code de programme pour executer les etapes d'un procede selon 
Tinvention, lorsque ledit programme est execute par un ordinateur. 

L'invention sera mieux comprise grace a la description qui va suivre donnee a 
titre d f exemple non limitatif en reference aux dessins annexes dans lesquels : 

25 - La figure 1 represente sch6matiquement une matrice sismique 

tridimensionrielle et une trace sismique locale continue. 

- La figure 2 represente sch6matiquement, un exemple de voisinage local N^k 
consistant en un ensemble de trace sismique locale continues S pq> k(t) voisines de la 
trace centrale de reference Syjc(t) elle mdme situee a la verticale du g6ophone Gy et 

30 interpolant les amplitudes sismiques autour de t=k. 
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- La figure 3 represente sch6matiquement, me section verticale de voisinage 
local conditionnel Nij,k(r) d f une trace sismique locale continue de reference Sy,k(t). 

- La figure 4 represente schematiquement, une section verticale d f une matrice 
sismique tridimensionnelle avec un decalage optimal h^k et un voisinage 

5 conditionnel N^r) . . 

- La figure 5 represente schematiquement, une section verticale d'une matrice 
sismique tridimensionnelle analogue a la figure 4 avec presence d'une faille ou 
discontinuity 

- La figure 6 represente schematiquement un organigramme fonctionnel d'un 
1 0 proc6de selon Tinvention. 

- La figure 7 represente schematiquement, une vue partielle agrandie d'une 
matrice sismique tridimensionnelle pr£sentant une discontinuite et illustrant la mise en 
oeuvre du procede selon l'invention. 

En reference aux figures 1 a 5, les elements identiques ou fonctionnellement 
15 equivalents sont designs ou reperes de maniere identique. 

Sur la figure 1, une matrice sismique tridimensionnelle est obtenue par releve 
de mesures enregistr6es par des g6ophones Gij disposes sur un r6seau x,y en des points 
de coordonnees i, j. Le releve des mesures echantillonnees dans le temps est repr6sente 
suivant un axe t descendant representatif de la profondeur ou d'une verticale 
20 descendante k partir de la surface de la terre ou de la mer. Les mesures sont 
caracterisees par leur amplitude, par exemple une amplitude relevee par le geophone 
Gy* au temps ou a la profondeur d'echantillonnage tk. La mesure discrete effectuee par 
le geophone Gij au temps ou k la profondeur t k est appelee amplitude sismique Sij >k . 

L'ensemble des amplitudes sismiques correspondant a un geophone Gij de 
25 coordonnees i, j est une matrice unidimensionnelle (S^, Sij 2 ,. . Sjjk,. - Syw) appeiees 
trace sismique discrete, car cette matrice unidimensionnelle correspond a la trace selon 
le point de coordonnees horizontals i, j de la matrice sismique tridimensionnelle 
obtenue par mesures sismiques. 

L'axe vertical t oriente selon une verticale descendante designe usuellement le 
30 temps, mais peut 6galement etre traite pour representer une profondeur a partir de la 
surface. 
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Uinvention concerne aussi bien Implication a une troisieme coordonnSe t 
representative du temps, qu!k une troisi&ne coordonn6e t representative de la 
profondeur. 

A partir de la trace sismique discrete situ6e k la verticale d'un geophone Gy on 
5 d&Bnit, par interpolation ou approximation des valeurs discretes autour de t=t k =k, une 
fonction continue Sij, k (t) qui est d6signee comme "trace sismique locale continue". Les 
m6thodes d f approximation ou ^interpolation de valeurs discretes pour engendrer une 
fonction continue sont nombreuses, et comprennent notamment les interpolations ou 
approximations polynomials, ainsi que les interpolations ou approximations 
1 0 trigonom&riques polynomials. 

Toute autre variante d'interpolation ou d'approximation fournissant une 
fonction continue peut egalement etre appliqu^e a la presente invention pour fournir 
une "trace sismique locale continue". 

Sur la figure 2, plusieurs traces sismiques locales continues definissent un 
15 voisinage d'une trace sismique locale continue "centrale" Sij, k (t) de reference. Le 
voisinage d f une trace sismique locale continue Sy,k(t) est defini comme Tensemble des 
traces sismiques locales continues dont les indices spatiaux horizontaux p,q sont 
voisins des indices spatiaux horizontaux i, j de la trace sismique locale continue de 
r6ference. 

20 A litre d ! exemple, les coordonnSes spatiales horizontales p,q correspondant au 

geophone G pq sont voisines des coordonn6es spatiales horizontales i, j correspondant 
au geophone Gy si les valeurs absolues des differences i-p et j-q sont inferieures k des 
entiers donn6s, par exemple k 2. 

Dans ce cas, comme represents k la figure 2, la trace sismique locale continue 
25 Sij,k(t) est associSe k huit traces sismiques locales continues voisines entourant la trace 
sismique locale continue "centrale" S ij>k (t). 

Dans le cas de traces sismiques locales continues produites par des mesures 
sismiques, le profil des horizons g6ologiques introduit des d6calages verticaux entre 
des traces sismiques locales continues voisines. Dans le but de d6terminer des relations 
30 entre deux traces sismiques continues locales voisines Si Jfk (t) et S pq>k (t) centres sur la 
meme coordonnSe verticale d'Schantillonnage t=t k =k et correspondant k des 
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coordonnees spatiales i, j et p,q difierentes, on calcule la fonction de correlation 
Rij^kCh) des deux traces sismiques locales continues. 

La fonction de correlation de deux traces sismiques locales continues voisines 
est obtenue par la formule suivante 

iW) — 



oil le num6rateur correspond a la fonction de covariance de Sij,k(t) et S pq ,k(t) 

tk+A 

obtenue par Texpression suivante Qj,pq,k(li) = |s ijk (t). S pq>k (t + h).dt . 

tk-A 

Dans cette integrate definissant C{j 9Vq ^(ti) , le parametre A d&finit une « fenetre verticale 
10 d'investigation » autour de t==t k ==k 

Par exemple, si les traces sismiques locales continues Sij,k(t) et S p(L k(t) sont des 
polyndmes trigonometriques de la forme suivante interpolant les donnes sismiques 

m 

S&kO) - £ a s M - cos ( s ^ 0 + b s ,k sin(sa> 1) 

s=l 
m 

S f j,k(t) = £a* k XOS(S^t) + b^sin(Stf>t) 



alors on peut montrer que la fonction de covariance Qj, pq ,k(h) est elle m6me un 
polynome trigonometrique de la forme suivante ou les coefficients A s lj,pq,k et B s 1J,pq>k 
dependent des coefficients a^* a.™*, b s ij,k et b s pq > k : 



20 C y , pq , k (h) = £ A* pq *\cos(stf> h) + Bi j ' pq ' k .sin(so> h) 



25 



Un resultat math£matique classique et connu est que la fonction de correlation 
Rii»ajc(li) = / XJ,pq,k traduit une similarite des fonctions correl6es lorsque 

cette fonction de correlation se rapproche de 1. 



L'etude des fonctions de correlation des traces sismiques locales continues 
Si,vk(t) et Spq,k(t) permet de definir comme decalage optimal la valeur hij, pq>k de h 
correspondant au maximum de la fonction de correlation Rij )Pq ,k(h) 3 c f est-&-dire au 



WO 2004/001450 



PCT/FR2003/001753 



7 



maximum de la fonction de correlation le plus proche de 1. Contrairement k Tart 
existant base sur une formulation discrete de la fonction de correlation R^pq^Oi), 
l'utilisation d'une formulation continue des traces sismiques et done de la fonction de 
correlation peimet d'obtenir un decalage optimum qui n'est pas contraint k etre un 
5 multiple entier du pas d'echantillonnage suivant I'axe vertical correspondant k la 
variable t 

La definition du decalage optimal hij,pq,k de deux traces sismiques locales 
continues voisines permet d'obtenir une premiere approximation des horizons passant 
par cette trace sismique locale continue de reference Sij,k(t), comme represent^ aux 
10 figures 3 k 5. Contrairement a Tart existant, le fait d'utiliser des decalages optimaux 
qui ne sont pas contraint k etre un multiple entier du pas d'echantillonnage suivant 
Paxe vertical permet d'6viter des erreurs connues sous le nom de «aliasing». 

Sur ces figures 3 a 5, on definit un voisinage conditionnel Nij,k(r) de la trace 
sismique continue locale de reference Sij >k (t) qui est un sous-voisinage du voisinage 
15 initial des traces sismiques locales continues retenues pour le calcul de correlation et 
de definition des decalages optimaux. Le voisinage conditionnel Nij,k(r) est choisi de 
maniere que, pour toute trace sismique continue locale Spq^t) appartenant k N ij>k (r), le 
decalage optimal hij >pq>k correspondant est tel que la correlation Rij^Chij^k) entre 
Sij,k(t) et S pq , k (t) est sup6rieure a un seuil predetermine r compris entre 0 et 1 . 

20 Comme on le voit sur la figure 5, la fixation du seuil r permet egalement 

d'englober dans un voisinage conditionnel une discontinue ou une faille geologique, 
ce qui constitue un avantage important par rapport k Tart anterieur. 

Les dispositions precitees fournissent ainsi ime modeiisation continue 
permettant la mise en ceuvre de Tinvention. 
25 En reference a la figure 6, un procede selon Tinvention comporte une 6tape de 

modeiisation numerique 100 pour d6finir des traces sismiques locales continues, des 
decalages optimaux et des voisinages conditionnels, comme decrit en reference aux 
figures 1 a 5. 



30 lignes et de colonnes correspondent aux coordonnees des geophones de mesures 
sismiques. 



A retape 101, on definit une matrice bidimensioimelle dont les indices de 



A retape 102, on choisit un point «graine» P(i j,t) ou premier point d'extraction. 
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A l'etape 103, on determine le noeud (ij,k) de la matrice sismique 
tridimensionnelle le plus proche du point «graine» P(i j,t). 

Ce noeud a pour coordonn6es spatiales i j et pour coordonnee temporelle ou de 
profondeur la coordonnee tk=k la plus voisine de la coordonn6e t du point «graine» 
5 P(i j,t) choisi k retape 102. 

A retape 104, on selectionne tout d'abord les traces sismiques continues locales 
S pq ,k(t) appartenant au voisinage conditionnel N^kCO defini k l'etape 100 et telles que 
les indices (p,q) correspondent a un emplacement vide de la matrice bidimentionnelle 
definie a I'etape 101. 

10 Ensuite, pour chaque trace sismique continue locale S pq> k(t) ainsi selectionnee, 

on attribue au point P(p,q,t) la valeur P(p,q,t + h ij>pq ,k) , ou h iifVqM est le d<§calage optimal 
du point P(p,q,t) par rapport au point voisin P(i j,k) determine a l'etape 103. Ce point 
P(P>q,t + hij, pq ,k) est consid6re comme un nouveau point situ6 sur Phorizon et est stocke 
dans Pemplacement (p,q) de la matrice bidimensionnelle defini k l'etape 101 . 

15 Le procede qui vient d'etre decrit est destine k determiner un ensemble de 

points appartenant a Phorizon passant par un point P(ij,t) , mais est egalement 
applicable pour attacher k ces points, au fur et a mesure de leur determination, un 
ensemble de propriet£s caract6risant la nature physique des terrains au voisinage de 
Phorizon ainsi determine. 

20 A cet effet, on peut prevoir d'attacher a chaque point P(p,q,t + hij, p ^k) de la 

matrice bidimensionnelle definie k Tetape 101 un ensemble de proprietes appeie 
"attributs sismiques". 

Ces proprietes sont gen6ralement representees par des couleurs differentes sur 
les visualisations d f ecran, de maniere k permettre une visualisation rapide et sans 
25 risque d'erreur des proprietes associees k un horizon geologique. 

En pratique, on effectue des etapes similaires aux etapes decrites en reference 
aux 6tapes 100 k 104 et Pon construit les attributs du point P(p,q,t + hij, pq> k) a Paide de 
Pequation de la trace sismique continue locale S pq>k (t) ; par exemple on construit ainsi 
un attribut enveloppe, un attribut phase ou un attribut frequence, de maniere connue en 
30 soi, par exemple d'aprfes Tenseignement du document de TANER M.T., KOEHLER 
F.,SHERIFF R.E., (1979), Complex seismic trace analysis, Geophysics, volume 44, 
n° 6. naeres 1 041 a 1063. 



WO 2004/001450 



9 



PCT/FR2003/001753 



Cette m6thode permet ainsi d'afficher Tattribut sismique sur l'horizon extrait en 
le "peignant" avec les couleurs correspondantes. 

Le proc6de venant d'etre d6crit permet ainsi de d6finir des valeurs d f extraction 
d'un horizon predetermine correspondant k une donnee r6elle, mais uniquement dans 
5 un voisinage conditionnel du point P(ij,k) le plus proche du point «graine» P(ij,t) 
choisi. 

Pour propager le proc£d6 d'extraction a Tensemble de la matrice 
bidimensionnelle d£finie k l'etape 101, on construit tout d'abord un ensemble de 
nouvelles graines constitue par les nouveaux points P(p,q,t + hij >pq) k) situe sur l'horizon 
10 qui viennent d'dtre determines et est stocks dans les emplacements (p,q) 
correspondants de la matrice bidimensionnelle defini a l'etape 101. Chacun de ses 
points est m6moris6 a l'6tape 106, pour Stre r6inject6 k l'etape 102, en y servant de 
nouveau point «graine» . 

Ceci propage le traitement successivement k Tensemble de la matrice 
15 bidimensionnelle correspondant aux geophones. 

Lorsque qu'aucune nouvelle graine ne peut etre trouv6e, le precede est oriente 
a l f 6tape 106 vers une etape 107 de visualisation sur 6cran des horizons ou des attributs 
ou propri6tes associ£es d'horizon g£ologique, de maniere k permettre une visualisation 
en couleur ou en niveaux de gris des proprietes associ6es a Thorizon geologique extrait 
20 correspondant au temps r6el t. 

Ainsi, a la fin de Talgorithme d6fini en reference k la figure 6, la matrice 
bidimensionnelle de l f 6tape 101, contient une plurality de points correspondant au 
meme horizon que le point initial «graine» choisit P(i j,t). 

L'ensemble des points peut 6tre reordonne avec ses voisins, conformement k la 
25 figure 7 pour definir un maillage recouvrant enti&rement l'horizon contenant le point 
«graine» P(i j,t), ce maillage permet ainsi de repr£senter l'horizon extrait comme une 
surface localement continue sur l'£cran de visualisation d'une station de travail. 

Le precede selon 1'invention est de preference mis en ceuvre sur un dispositif 
comportant des moyens appropries a la mise en oeuvre des etapes decrites en reference 
30 a la figure 6. 
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En particulier, un dispositif selon l'invention comporte des moyens de 
memorisation n6cessaires aux calculs successifs et les moyens de visualisation 
n6cessaires aTetape 107. 

De preference, un dispositif selon Tinvention est un dispositif programmable 
5 commande par un produit-programme d'ordinateur comportant des elements de code 
de programme pour ex6cuter les Stapes du proc6d6 d6crit en reference k la figure 6. 

L'invention d6crite en r&f&rence a plusieurs objets particuliers n'y est nullement 
limitee, mais couvre au contraire toute modification de forme et toute variante de 
realisation dans le cadre et l f esprit de l'invention. 
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REVENDICATIONS 

1. Proced6 pour extraire un horizon g£ologique et des propri6t6s associ6es d'une 
representation sismique, dans lequel on construit vine fonction continue Sij,k(t) par 
interpolation ou approximation des traces sismiques discretes d'une matrice sismique 

5 multidimensionnelle, ladite fonction 6tant d6sign6e comme "trace sismique locale 
continue", comportant les Stapes suivantes : 

a) , utiliser comme d6calage (vertical) optimal de deux traces sismiques locales 
continues voisines, la valeur de decalage rendant maximale leur fonction de 
correlation, ce decalage optimal n'&ant pas contraint k etre un multiple entier du pas 

10 d'6chantillonnage vertical ; 

b) retenir comme voisinage conditionnel d'une trace sismique locale continue 
"centraie" Sij >k (t) le sous-voisinage consistant en des traces voisines S pq>k (t) 
correspondant a des decalages optimaux hij >pq>k associes a des correlations Rij, pq ,k(hij,pq,k) 
sup6rieures h un seuil predetermine compris entre 0 et 1 ; 

15 c). constraire une matrice bidimensionnelle d'extraction destinee a §tre remplie 

par des points extraits appartenant au meme horizon que celui passant par le point 
graine ; 

d). choisir un point «graine» P(i j,t) et determiner le point P(ij,k). de la matrice 
sismique tridimensionnelle verticalement le plus proche ; 

20 e). estimer les propri6tes associees du voisinage conditionnel et remplir la 

matrice bidimensionnelle d'extraction par des proprietes d£cal6es par translation de la 
variable courante (t) de la valeur du decalage optimal (hij fP q^) correspondant au point 
(i, j, k) verticalement le plus proche. 

2. Precede selon la revendication 1, dans lequel on utilise comme point 
25 «graine» de Tetape d) tous les points nouvellement stock£s dans la matrice 

bidimensionnelle de Tetape b) et non encore utilises comme point graine. 

3. Proc6d6 selon la revendication 1 ou 2, dans lequel on remplit successivement 
le contenu de la matrice bidimensionnelle de l'etape c) au cours des iterations 
successives d'extraction. 
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4. Precede selon la revendication 1 ou 2, dans lequel on remplace 
successivement le contenu de la matrice bidimensionnelle de l f 6tape c) par une 
moyenne des contenus successifs au cours des iterations successives d'extraction. 

5. Precede selon la revendication 1 ou 2, dans lequel la propriety de subsurface 
5 extraite est un ensemble d' attributs sismiques calculus en chaque point extrait sur 

l'horizon passant par le point graine, le calcul de ces attributs etant lui meme realise 
simultanement k V extraction de ces points. 

6. Precede selon la revendication 1 ou la revendication 5, dans lequel on 
visualise sur un ecran de visualisation les attributs sismiques « peints » sur Phorizon 

10 extrait. 

7. Dispositif pour la mise en ceuvre du precede selon Tune quelconque des 
revendications 1 k4 y comportant des moyens pour utiliser comme decalage optimal de 
deux traces sismiques locales continues voisines, la valeur de decalage rendant 
maximale leur fonction de correlation, des moyens pour retenir comme voisinage 

15 conditionnel d'une trace sismique locale continue "centrale" de reference Sij,k(t) le 
sous-voisinage consistant en des traces voisines S pqu k(t) correspondant k des decalages 
optimaux hij,pq,k associes a des correlations Rij,pq^(hij, P q^) superieures k un seuil 
predetermine compris entre 0 et 1, des moyens pour construire une matrice 
bidimensionnelle d'extraction destinee a 8tre remplie par des valeurs extraites, des 

20 moyens pour choisir un point «graine» P(ij,t) et determiner le point P(ij,k) 
verticalement le plus proche et des moyens pour estimer les propriet6s associees du 
voisinage conditionnel et remplir la matrice bidimensionnelle d'extraction par des 
proprietes decaiees par translation de la variable courante (t) de la valeur du decalage 
optimal (hy^k) correspondant au point (i, j, k) verticalement le plus proche. 

25 8. Dispositif selon la revendication 7, comportant des moyens de memorisation 

et des moyens de visualisation de parametres sismiques determines k Taide du precede 
selon Tune quelconque des revendications 1 a 6. 

9. Produit-programme d'ordinateur, comportant des elements de code de 
programme pour executer les etapes du precede selon Tune quelconque des 
30 revendications 1 & 5, lorsque ledit programme est execute par un ordinateur. 
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10. Produit-programme d'ordinateur, comportant des elements de code de 
programme pour exScuter les Stapes du procede selon la revendication 6, lorsque ledit 
programme est execute par un ordinateur. 
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Dale de priori te (jour/mois/annee) 
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1. 11 est notifie par la presenle qu'a la date indiquec ci-dessus comme date d'expedition de cet avis, le Bureau international a communique, 
comnie le prevoit Particle 20, la demande intemationale aux offices designes suivants : 

AU, A2, BY, CH, CN, CO, DE, DZ, EP, HU, JP, KG, KP, KR, MD, MK, MZ, RU, TM, US 

Confomiement a la regie 47.1.C), troisieme phrase, ccs offices acceptent le present avis comme preuve determinante du fait que la 
communication de la demande inlernationale a bion eu lieu a la date d'expedition indiquee plus haul, el le deposani n'est pas tenu de 
remettre de copie de la demande inlernationale a I' office ou aux offices designes. 

2. Les offices designes suivants ont renonc6 a 1' exigence selon laquelle cette communication doit etre effectuee a cette dale: 

AE, AG, AL, AM, AP, AT, BA, BB, BG, BR, BZ, CA, CR, CU, CZ, DK, DM, EA, EC, EE, ES, Fl, GB, GD, GE, GH, GM, 
HR, ID, IL, IN, IS, KE, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MG, MN, MW, MX, Nl, NO, NZ, OA, OM, PH, PL, PT, RO, 
SC, SD, SE, SG, SK, SL, TJ, TN, TR, TT, TZ, UA, UG, UZ, VC, VN, YU, ZA, ZM, ZW 

La communication sera effectuee seulement sur demande de ces offices. De plus, le deposant n'est pas tenu de remettre de copie de la 
demande inlernationale aux offices en question (regie 49.l)a-bis)). 

3. Le present avis est accompagne d'une copie de la demande intemationale publiee par le Bureau international le 
31 decembre 2003 (31.12.2003) sous le numero WO 20047001450 

4. DELAIS pour Ja presentation d'une demande d'examen preliminaire international et pour 1'ouverture de la phase nationale 

Le delai applicable pour 1'ouverture de la phase nationale sera, sous reserve de ce qui est dit au paragraphe suivant, de 30 MOIS a 
compter de la dale de priority, non seulement en ce qui conceme tout office elu lorsqu'une demande d'examen preliminaire international 
aura ete presentee avant l'expiralion du delai de 19 inois a compter de la date de priorite (voir 1'article 39.1)), mais egalement en ce qui 
concerne tout office designe, en r absence de presentation d'une telle demande d'examen, lorsque 1'article 22.1) lei que modifie avec effet 
au lcr avril 2002 sera applicable audit office designe. Pour plus de renseignements, voir la Gazette du PCT no 44/2001 du ler novembre 
2001, pages 19927, 19933 et 19935, ainsi que le bulletin PCT Newsletter, numeros d'octobre el de novembre 2001 et de fevrier 2002. 

En pratique, des delais autres que celui de 30 inois vont continuer de s'appliquer, pour des durees variables, en ce qui concerne certains 
offices designes el elus^ Pour des mises a jour regulieres quant aux delais applicable^ (20, 2 1, 30 ou 31 mois ou autre delai), office par 
office, on se reportera a la Gazette du PCT, au bulletin PCT Newsletter ainsi qu'aux chapitres nationaux pertinents dans le volume II du 
Guide du deposant du PCT, accessibles sur le site Internet de I'OMPl, par 1' intermedia ire de liens a partir de diverses pages du site, y 
compris celles de la Gazette, de. la Newsletter et du Guide, a 1'adresse suivante : http://www.wipo.int/pcL/fr/index.html. 

Quant a la presentation d'une demande d'examen preliminaire international, voir le Guide du deposant du PCT, volume I/A, chap i Ire 
IX. Seul un deposant qui est ressortissanl d'un Etat contractanl du PCT lie par le chapitre II ou qui y a son domicile peut presenter unc 
demande d'examen preliminaire international (actuel lenient, tous les Etats coniractants du PCT sont lies par le chapitre 11). 

Le deposant est seul responsable du respect de tous les delais vises ci-dessus. 
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| X | le deposant X I'inventeur | | le mandataire j ~[ le representant commun 


Nom et adresse 
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Nationality (nom de I'Etat) 
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2. Le Bureau international notifie au d&posant que le changement indique ci-apres a 6te enregistre en ce qui concerne: 
| | la personne | I le nom j | Tadresse [ ^ la nationality [ ] le domicile 


Nom et adresse 

BOSQUET, Fabien 
33, rue de la Ravinelle 
F-54000 NANCY 
FRANCE 


Nationality (nom de I'Etat) 
FR 


Domicile (nom de I'Etat) 
FR 


no de telephone 


no de telecopieur 


no de telelmprimeur 


3. Observations complementaires, le cas Scheant: 


4. Une copie de cette notification a et6 envoyee: 
| X | a I'office recepteur j ^ aux offices designes concerned 
j | a I'administration chargee de la recherche internationale "j" J aux offices elus concerned 
[ | a I'administration chargee de I'examen preliminaire international | | autre destinataire: 
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(PCT Article 36 and Rule 70) 


Applicant's or agent's file reference 
T020575 JKLC 


-» TOrrm?T> . __ T __ T See Notification of Transmittal of International 
FOR FURTHER ACTION Pre iiminary Examination Report (Form PCT/IPEA/416) 


International application No. 


International filing date (day/month/year) 


Priority date (day/month/year) 


PCT/FR2003/001753 
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19juin 2002(19.06.2002) 


International Patent Classification (IPC) or national classification and IPC 
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This international prelirninary examination report has been prepared by this International Preliminary Examining Authority 
and is transmitted to the applicant according to Article 36. 



This REPORT consists of a total of 



. sheets, including this cover sheet. 



I 1 This report is also accompanied by ANNEXES, i.e., sheets of the description, claims and/or drawings which have been 
' — ' amended and are the basis for this report and/or sheets containing rectifications made before this Authority (see Rule 
70. 1 6 and Section 607 of the Adrninistrative Instructions under the PCT). 



These annexes consist of a total of 



sheets. 



3. This report contains indications relating to the following items: 
Basis of the report 
Priority 

Non-establishment of opinion with regard to novelty, inventive step and industrial applicability 
Lack of unity of invention 

Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

Certain documents cited 

Certain defects in the international application 

Certain observations on the international application 



I 


I2SI 


II 


□ 


m 


□ 


IV 


□ 


V 




VI 


□ 


VII 


□ 


vra 


□ 



Date of submission of the demand 

24 novembre 2003 (24.1 1 .2003) 


Date of completion of this report 

23 March 2004 (23.03.2004) 


Name and mailing address of the IPEA/EP 
Facsimile No. 


Authorized officer 
Telephone No. 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



Im 



al application No. 

PCT/FR2003/001753 



I. Basis of the report 



1. With regard to the elements of the international application:* 
the international application as originally filed 
the description: 

pages 

pages 



1-10 



, as originally filed 
, filed with the demand 



, filed with the letter of 



the claims: 



1-10 



pages 



, as originally filed 

, as amended (together with any statement under Article 19 
, filed with the demand 



filed with the letter of 



the drawings: 

pages 

pages 

pages 



1/4-4/4 



_ , as originally filed 
, filed with the demand 



I I the sequence listing part of the description: 



mt filed with the letter of 



pages 



, as originally filed 

_ , filed with the demand 



filed with the letter of 



2. With regard to the language, all the elements marked ahove were available or furnished to this Authority in the language in which 
the international application was filed, unless otherwise indicated under this item. " 

These elements were available or furnished to this Authority in the following language which is* 

□ 

the language of a translation furnished for the purposes of international search (under Rule 23.1(b)). 

□ the language of publication of the international application (under Rule 48.3(b)). 

the language of the translation furnished for the purposes of international preliminary examination (under Rule 55.2 and/ 
or 55.3). 

3. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international 
preliminary examination was carried out on the basis of the sequence listing- 

□ contained in the international application in written form 

□ filed together with the international application in computer readable form. 

□ furnished subsequently to this Authority in written form 

□ furnished subsequently to this Authority in computer readable form. 

The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished. 

The statement that the information recorded in computer readable form is identical to the written sequence listing has 
been furnished. 



□ 



The amendments have resulted in the cancellation of: 

□ the description, pages 

the claims, Nos. 

the drawings, sheets/fig 



5. n ™ S rep ° rt haS beCn estabnslied 33 if (some of) the amendments had not been made, since they have been considered to go 
beyond the disclosure as filed, as indicated in the Supplemental Box (Rule 70.2(c)).** 

* Replacement sheets which have been furnished to the receiving Office in response to an invitation under Article 14 are referred to 
m rli r . e J? ort 05 " ori Sinally filed" and are not annexed to this report since they do not contain amendments (Rule 70.16 
and 70 J 7). * y 

**Any replacement sheet containing such amendments must be referred to under item 1 and annexed to this report. 
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V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 



I. Statement 

Novelty (N) 

Inventive step (IS) 
Industrial applicability (IA) 



Claims 
Claims 

Claims 
Claims 

Claims 
Claims 



1-10 



1-10 



1-10 



YES 
NO 
YES 
NO 

YES 
NO 



Citations and explanations 

1. Reference is made to the following documents ; 

Dl: US-A-5 615 171 

D2: US-A-6 138 076 

D3: US-A-4 672 546 

D4: US-A-5 930 730 



1.2 D5 (EP 0447249) was not cited in the international 
search report. A copy of said document is attached. 

2. Claim 1 discloses a method for extracting a 
geological horizon; whereas claims 7, 9 and 10 
disclose , respectively, a device and program 

products for performing the steps of said method of 

>» 

claim 1. 

3. Dl, which is considered to be the prior art closest 
to the subject matter of claim 1, describes a method 
for extracting seismic horizons .that uses a seed 
point in a seismic trace and subsequently performs a 
comparison or a discrete analysis of the correlation 
between the adjacent traces to establish geological 
horizons . 
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PCT 




application No. 
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5. 



5.1 



6. 



With regard to the closest prior art, the subject 
matter of claim 1 differs from Dl in that a 
continuous function is constructed by interpolating 
or approximating discrete seismic traces and an 
optimum offset value, which is not necessarily an 
integer multiple of the vertical sampling pitch, is 
used. The seismic traces are compared by determining 
the optimum offset value associated with 
correlations greater than a predetermined threshold. 
The subject matter of claim 1 is therefore novel 
(PCT Article 33 (2) ) . 

The problem that the present invention is intended 
to solve can be considered to be that of reducing or 
eliminating the aliasing effect. 

The solution proposed in claim 1 of the present 
application is considered inventive (PCT Article 
33(3)), because the combination of the features set 
forth in claim 1 is not a routine technical measure 
for a person skilled in the art and none of 
documents D2 to D5 suggests the solution proposed in 
the present application. 

Claims 2 to 6 and 8 are dependent, respectively , on 
claims 1 and 7 and describe several ways and options 
available for implementing the invention according 
to claim 1. Said dependent claims therefore also 
comply, as such, with the PCT requirements of 
novelty and inventive step. 

The industrial applicability of the method and the 
device described in claims 1 and 7 is obvious. 
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(article 36 et regie 70 du PCT) 



Reference du dossier du deposant ou du 
mandataire 


POUR SUITE A DONNER vo £ ,a notification de transmission du rapport d'examen 

preliminaire International (formulaire PCTVIPEA/416) 


Demande internationale No. 
PCT/FR 03/01753 


Date du depdt international Qour/mols/ann6e) 
11.06.2003 


Date de priority QourAnoIsfann6e) 
19.06.2002 


uiassmcatton Internationale des brevets (CIB) ou a la fois classification Rationale et CIB 

G01V1/30 


Deposant " " . 

EARTH DECISION SCIENCES 



1. 



2. Ce RAPPORT comprend 5 feuilles, y compris la presente feuille de couverture. 



1=1 e? 'SSSZSSi^ d ! ,eU i" eS dS ' a desc "P« on . des revendications ou des dessins qui 

Ces annexes comprennent feuilles. 



Le present rapport contient des indications et les pages correspondantes relatives 



aux points suivants : 



1 




II 


□ 


III 


□ 


IV 


□ 


V 




VI 


□ 


VII 


□ 


VIII 


□ 



Base de I'opinion 
Priority 

JSS^ qUant * !a nouve aute, I'activite inventive et la 

possibility duplication industrielle 

Absence d'unite de invention 

d!2^ 66. 2 (a)(ii) quant a la nouveaute, I'activite inventive et la possible 

d application industrielle; citations et explications a I'appul de cette declaration 

Certains documents cites 

Irregularites dans la demande Internationale 

Observations relatives a la demande Internationale 



Date de presentation de la demande d'examen preliminaire 
internationale 



24.1 1 .2003 



Norn et adresse postale de I'admlnstration charges de I'examen 
preliminaire international 



Office europeen des brevets - P.B. 5818 Patenttaan ? 
NL-2280 HV Rijswljk - Pays Bas «nuaan 2 

Tel. +31 70 340 - 2040 Tx: 31 651 epo nl 
Fax:+31 70 340-3016 



Date d'achevement du present rapport 
23.03.2004 



Fonctionnalre autorise 
Modesto, C 

N° de telephone +31 70 340-1055 




= ormulaire PC771PEA/409 (feuille de titre) (Janvier 2004) 



RAPPORT D'EXAMEN 
PRELIMINAIRE INTERNATIONAL 



Demande Internationale n° PCTFR 03A)1 753 



I. Base du rapport 

1 . En ce qui concerne les elements de la demande Internationale (les feuilles de remplacement qui ont et6 
remises a Voffice recepteur en reponse a une invitation faite conformement it I'article 14 sont consider6es, dans 
le present rapport , comme "initialement deposees" et ne sont pas jointes en annexe au rapport puisqu'elles ne 
contiennent pas de modifications (regies 70. 16 et 70. 17)) : 

Description, Pages 

1-10 telles qu'initialement deposees 
Revendications, No. 

1 -1 0 telles qu'initialement deposees 
Dessins, Feuilles 

1/4-4/4 telles qu'initialement deposees 

2. En ce qui concerne la langue, tous les elements indiques ci-dessus 6taient a la disposition de I'administration 
ou lui ont ete remis dans la langue dans laquelle la demande internationale a ete deposee, sauf indication 
contraire donn§e sous ce point. 

Ces elements 6taient k la disposition de I'administration ou lui ont et6 remis dans la langue suivante: ,qui est: 

□ la langue d'une traduction remise aux fins de la recherche internationale (selon la regie 23.1(b)). 

□ la langue de publication de la demande internationale (selon la rfegle 48.3(b)). 

□ la langue de la traduction remise aux fins de I'examen preliminaire internationale (selon la reqle 55.2 ou 
55.3). 

3. En ce qui concerne les sequences de nucleotides ou d'acide aminesdivulgu6es dans la demande 
internationale (le cas 6cheant), I'examen preliminaire internationale a «§t§ effectue sur la base du listage des 
sequences : 

□ contenu dans la demande internationale, sous forme ecrite. 

□ d§pos6 avec la demande internationale, sous forme d§chiffrable par ordinateur. 

□ remis ulterieurement k I'administration, sous forme ecrite. 

□ remis ulterieurement k I'administration, sous forme dechiffrable par ordinateur. 

□ La declaration, selon laquelle le listage des sequences par 6crit et fourni ulterieurement ne va pas au-dela 
de la divulgation faite dans la demande telle que d6pos6e, a et§ fournie. 

□ La declaration, selon laquelle les informations enregistn§es sous dechiffrable par ordinateur sont identiques 
k celles du listages des sequences Presente par ecrit, a et6 fournie. 

4. Les modifications ont entraTn6 I'annulation : 

□ de la description, pages : 

□ des revendications, nos : 

□ des dessins, feuilles : 
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5 ' D tl^T 0 !! tff f |? rmUl6 ^straction faite (de certaines) des modifications, qui ont ete considerees 

70 2(c)) • XP ° Se ,nvention tel ^ a 616 d6 P° s6 - comm e ■ est indique ci-apres (regie 

etann^ dGS modifications de cette nature d °'t etre indiquee au point 1 

6. Observations complementaires, le cas echeant : 



V. Declaration motivee selon I'article 35(2) quant a la nouveaute, I'activite inventive et la possibility 
d'application industnelle; citations et explications a I'appui de cette declaration 

1. Declaration 
Nouveaute 

Activite inventive 



Oui: 


Revendications 


1-10 


Non: 


Revendications 




Oui: 


Revendications 


1-10 


Non: 


Revendications 




Oui: 


Revendications 


1-10 


Non: 


Revendications 





2. Citations et explications 
voir feuille separee 
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Concernant le point V 

Declaration motivee selon I'article 35(2) quant a la nouveaute, I'activite inventive et 
la possibility d'application industrielle; citations et explications a I'appui de cette 
declaration. 

1 II est fait reference aux documents suivants: 

D1: US-A-5 615 171 
D2: US-A-6 138 076 
D3: US-A-4 672 546 
D4: US-A-5 930 730 

1.2 Le document D5: EP 0447249 n'a pas ete cite dans le rapport de recherche 
international. Une copie de ce document est joints en annexe. 

2 La revendication 1 devoile un procede pour extraire un horizon geologique, les 
revendications 7, 9 et 10 devoilant respectivement un dlspositif et des produits- 
programmes pour executer les etapes du dit procede de la revendication 1 . 

3 Le document D1 , qui est considere comme etant I'etat de la technique le plus proche 
de I'objet de la revendication 1, decrit un procede pour extraire des horizons 
sismiques lequel utilise un point grain dans une trace sismique et apres execute une 
comparaison ou une etude discrete de la correlation entre des traces sismiques 
voisines pour etablir des horizons geologiques. 

4 Vis a vis de I'etat de la technique le plus proche, I'objet de la revendication 1 differe 
de D1 en ce que: une fonction continue est construite par interpolation ou 
approximation des traces sismiques discretes et une valeur de decalage optimal qui 
n'est pas contrainte a etre un multiple entier du pas d'echantillonnage vertical est 
utilisee. La comparaison entre des traces sismiques est realise en determinant 
trouvant la valeur de decalage optimum associee a des correlations superieures a 
un seuil predetermine. L'objet de la revendication 1 est done nouveau (article 33(2) 
PCT). 



Le probleme que se propose de resoudre la presente demande peut done etre 
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considere comme 6tant la reduction ou I'elimination de I'effet "d'aliasing". 

5.1 La solution proposee dans la revendication 1 de la presente demande est consideree 
comme inventive (article 33(3) PCT) parce que I'ensemble des caracteristiques 
exposees dans la revendication 1 ne releve pas d'une demarche technique normale 
pour la personne du metier et, aucun des documents (D2 - D5) suggere la solution 
presente par la presente demande. 

6 Les revendications 2-6 et 8 dependantes respectivement des revendications 1 et 7 
decrivent plusieurs manieres et options pour realiser I'invention selon la revendication 
1 et satisfont done egalement, en tant que telles, aux conditions requises par le PCT 
en ce qui conceme la nouveaute et I'activite inventive. 

7 Lapplication industrielle du procede et du dispositif decrite dans les revendications 
1 et 7 est consideree evidente. 
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